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Abstract. We consider supersymmetric model with right-handed (s)neutrinos where the neutrino 
masses are purely Dirac-type. We discuss cosmology based on such a scenario, paying particular 
attention to the case that the gravitino is the lightest superparticles (LSP) while the right-handed 
sneutrino is the next-LSP. It will be shown that the cosmological constraints on the gravitino-LSP 
scenario (in particular, those from the big-bang nucleosynthesis) are drastically relaxed in such a 
case. We will also consider the implication of such scenario to the structure formation. 

PACS. 12.60.Jv Supersymmetric models - 95.35.-fd Dark matter - 98.80.Cq Particle-theory and 
field-theory models of the early Universe - 14.60.Pq Neutrino mass and mixing 



1 Introduction 

Existence of dark matter (DM) in our universe, which 
is strongly supported by a lot of recent cosmological 
observations [Tl|2l[3] , requires physics beyond the stan- 
dard model (SM). Many possibilities of dark matter 
have been discussed in various frameworks of particle 
physics models so far 4 . Importantly, properties of the 
dark matter particle depend strongly on the particle 
physics model we consider. 

In the framework of supersymmetry (SUSY) , prob- 
ably most popular candidate of dark matter is ther- 
mally produced lightest neutralino which is usually 
assumed to be the lightest superparticle (LSP). How- 
ever, if we try to build a supersymmetric model which 
accommodates with all theoretical and experimental 
requirements, we expect that there exist new exotic 
would-be-DM particles which are not superpartners 
of the SM particles, for example, gravitino ip^ and 
right-handed sneutrino which are superpartners of 
graviton and right-handed neutrino, respectively. If the 
neutrino masses are Dirac-type, in particular, v^s are 
expected to be as light as superpartners of SM parti- 
cles in the framework of the gravity-mediated SUSY 
breaking. Existence of these exotic superparticles may 
significantly change the phenomenology of dark matter 
in supersymmetric models. 

In this paper, we consider the supersymmetric model 
in which the neutrino masses are Dirac-type and dis- 
cuss cosmological implications of such a scenario. It 
has already been pointed out that the the i^i^-DM sce- 
nario can be realized in Ref. Here, we consider an- 
other case where the gravitino is the LSP and one of 
the right-handed sneutrinos is the next-to-the-lightest 
superparticle (NLSP) [6]. If the gravitino is the LSP, 



it may be a viable candidate of dark matter also in 
the case without the right-handed sneutrinos 7,8,9]. 
In such a case, however, stringent constraints on the 
scenario are obtained from the study of the gravitino 
production at the time of the reheating after inflation 
and also from the study of the big-bang nucleosynthe- 
sis (BBN) reactions. With the z^ij-NLSP, we reconsider 
cosmological constraints on the "i/j^-LSP scenario. We 
pay particular attention to the BBN constraints and 
also to the constraints from the structure formation of 
the universe. We will see that the BBN constraints are 
significantly relaxed if there exists the z/^-NLSP. 



2 Model Framework 

In this section, we discuss the model in which gravitino 
is the LSP while right-handed sneutrino is the NLSP. 
We assume that neutrino masses are purely Dirac- 
type, and the superpotential of the model is written 
as 

W = VTmssm + y^LHuO'j^, (1) 

where Wmssm is the superpotential of the minimal su- 
persymmetric standard model (MSSM), L = {OltBl) 
and Hu = (iJ+,iJ°) are left-handed lepton doublet 
and up- type Higgs doublet, respectively. In this ar- 
ticle, "hat" is used for superfields, while "tilde" is for 
superpartners. Generation indices are omitted for sim- 
plicity. 

In this model, neutrinos acquire their masses only 
through Yukawa interactions as rrii^ — y^{H'^) — y^v sin 
where v is the vacuum expectation value (VEV) of the 
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standard model Higgs field (i; ~ 174 GeV) and tan/3 = 
{H^) / {H^). Thus, the neutrino Yukawa coupling is de- 
termined by the neutrino mass through the equation: 
sin/3 = 3.0 x IQ-^^ ^ (m2/2.8 x IQ-^ eY^f/^. Mass 
squared differences among neutrinos have already been 
determined accurately at neutrino oscillation experi- 
ments p^miTT]. In this article, we assume that the spec- 
trum of neutrino masses is hierarchical, hence the largest 
neutrino Yukawa coupling is of the order of 10~^^. 
We use yi, = 3.0 x 10^^^ for the numerical analysis 
in this article. With right-handed (s)neutrinos, new 
soft SUSY breaking terms are introduced in addition 
to the usual terms of the MSSM. Those are right- 
handed sneutrino mass terms and tri-linear coupling 
terms called Ai^-terms. Breaking terms relevant to our 
analysis are 

-CsoFT = -Af|Z^L - M?^i>'^i>ii 

+ iA,LH^i^'i^ + h.c), (2) 

where all breaking parameters, M|^, Mj^j^ and A^, are 
defined at the electroweak (EW) scale. We parameter- 
ize Ai, by using the dimensionless constant a^, as A^ = 
a^y^Mj^. We adopt gravity-mediated SUSY breaking 
scenario and, in such a case, is expected to be 0(1). 
Though the A^-term induces the left-right mixing in 
the sneutrino mass matrix, the mixing is safely ne- 
glected in the calculation of mass eigenvalues due to 
the smallness of neutrino Yukawa coupling constants. 
Thus, the masses of sneutrinos are simply given by 
m?^ = M| -I- i cos(2/3)m|, m?^ = M?^, where mz ~ 
91 GeV is the Z boson mass. In the following discus- 
sion, we assume that all the right-handed sneutrinos 
are degenerate in mass for simplicity. 

In this article, we consider the V'm"LSP scenario 
with i'fl^-NLSP. In such a case, the next-to-next-LSP 
(NNLSP) plays an important role in the thermal his- 
tory of the universe. However, there are many possibil- 
ities of the NNLSP, depending on the detail of SUSY 
breaking scenario. Thus, we concentrate on the case 
that the NNLSP is the lightest neutralino, whose com- 
position is Bino B. This situation is easily obtained if 
we consider the so-called constrained-MSSM type sce- 
nario. It is not difficult to extend our discussion to the 
scenario with other NNLSP candidate. 



3 Constraints from BBN 

It is well known that models with the i/j^-LSP usually 
receive stringent constraints from the BBN scenario. In 
these models, LSP in the MSSM sector (which we call 
MSSM-LSP) is long-lived but decays to gravitino with 
hadrons, which may spoil the success of BBN scenario. 
The BBN constraints give the upper bound on Yx^^vis 
as a function of tx , where X stands for a long-lived but 
unstable particle, Yx = [nx / s]t<^Tx (with nx and s 
being the number density of X and the entropy density 
of the universe, respectively), E^is is the mean energy 
of visible particles emitted in the X decay, and tx is 
the lifetime of the particle X [12]. We use the upper 



bound on Yx-E-vis obtained in the studies. In order to 
evaluate BBN constraints quantatively, we calculate 
Bha^dYxEvis as a function of r^, then search for allowed 
paramter region. 

First, we wiU see the decay of the MSSM-LSP, 
which is i?-like neutralino, in order to calculate Shad^'vif 
and . Main decay modes are the following two-body 
decays: B — > vfiPl, B — > V'/i7> and B — > if^^Z. By the 
use of decay widths of these processes in Ref. [l^ . Im- 
portanly, the decay mode B — > vrvl competes with 
the mode B — > ■0^7 or it even dominates in total decay 
in wide parameter region, especially when the grav- 
itino mass 7713/2 is larger than 0.1 GeV when mofi = 
100 GeV and a„ = 1. We have also checked that the 
lifetime of B is 10^~^ seconds in that parameter region. 

Without the 7>fl-NLSP, the decay mode B ipul/Z 
dominates in total decay, and many visible particles 
are emitted through photon or Z boson, which may 
spoil the success of the BBN scenario. As a result, the 
gravitino mass is strictly constrained as 7773/2 ^0.1 
GeV for < 1 second [13]. In our scenario with the 
£'fl-NLSP, however, less hadrons are emitted, though 
the Bino-like neutralino is long-livedEl. Therefore, con- 
straints from BBN is expected to be relaxed in the 
"^3/2 ^ 0-1 GeV region. 

As we saw that no hadrons are emitted in two-body 
decays, in order to calculate i?had£'vis2 we consider 
three- or four-body decays: B i'fj.ff, B — > VRelfJ', 
and B — > i^Ri^Lff, where / and / denote fermion and 
anti-fermion, respectively. Although branching ratios 
of these processes are much smaller than 1, they have 
impacts on the BBN scenario. 

Lasty, we determine Yg by the use of following 
formula for (would-be) density parameter of B [13] : 
J7^/i^ = Cmodei X 0.1 [7)7-^/100 GcV] , whcrc 777^ is 
Bino mass and the additional parameter Cmodei is in- 
troduced to take the model dependence into account in 
our analysis: Cmodei ~ 1 for the neutralino in the bulk 
region, Cnodei '^0.1 for that in the co-annihilation or 
funnel region, and Cmodei ^10 for the pure Bino case 
without co-annihilation. 

Our numerical results are shown in Fig[Tl where 
the BBN constraints are depicted on the (7^3/2, tr^) 
plane. We take Cmodei — 10, 1,0.1 in the left, middle, 
and right figures, respectively. Other parameters are 
set as 77ip^ = 100 GeV, Oi, = 1, and 777;;^ = 1.5mj^. 
Shaded regions are ruled out by the BBN scenario. As 
shown in these figures, the constraints are drastically 
relaxed compared to those in models without the t/^- 
NLSP. 

As shown in the figures, new allowed region ap- 
pears; for example, for Cmodei = 1, 0.1 GeV < 7)7.3/2 < 
40 GeV. In that region, the BBN constraints give the 
upper bound onmg. Since the decay mode B ii^qq 
is subdominant in this region, this bound comes mainly 



Left-handed neutrinos injected by the decay might pos- 
sibly change the abundance of *He [T?]. However, we have 
checked that the BBN constraints on the neutrino injection 
are much less stringent than those on three- or four-body 
decays as shown in the following discussion. 
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Fig. 1. Constraints from tlie BBN on the (7713/2, rrig) plane: Parameters are cliosen to be mp^ — 100 GeV, a,^ = 1, 
and mfj^ = 1.5mg. We set Cmodoi =10, 1, and 0.1 in tiie left, middle, and right figures, respectively. The middle shaded 
regions are ruled out by the BBN scenario, while dark and light shaded regions are excluded by the WMAP measurement 
and the structure formation of the universe, respectively. 



from four-body decays, B vnej^qq' and B VRU^qq. 
This fact can be understood intuitively: Shad and -Evis, 
are enhanced when is large. On the contrary, in 
the 0.01 GeV < mg/a < 0.1 GeV region, Bino-like 
neutralino decays mainly into the gravitino through 
the B — > ■!/;^7 process with the lifetime Tj^ <\ second. 
Since the decay occurs before the BBN starts, it does 
not affect the BBN scenario. This situation also holds 
in the usual -(/;^-LSP scenario without the I'ij-NLSP, 
and the same allowed region can be seen in Ref.|13|. 
In the case of Cmodei = 10(0.1) in the left (right) fig- 
ure, the constraints from the BBN scenario is more 
(less) stringent than the Cmodci — 1 case. As a result, 
the upper bound on nig becomes smaller (larger) . We 
also find that the region rrig ~ 100 GeV with 0.1 GeV 
^ "^3/2 ^ 1 GeV is excluded in the left and mid- 
dle figures: this is because the process B vrvl is 
kinematically suppressed and branching ratio of the 
process B ip^j^Z is relatively enhanced. 

In addition to the BBN constraints, we also de- 
pict other cosmological bounds in Fig[l] the gravitino 
abundance originating in B must not exceed the value 
observed in the WMAP, /2dm/i^ ^ 0.105 [3], which is 
shown as a dark shaded region in Fig[T] This constraint 
gives the upper bound on 7713/2- Another constraint, 
•^3/2 0-4^DM, is also depicted as a light shaded re- 
gion, which comes from the structure formation of the 
universe as discussed in the next section. 



4 Constraints from Structure Formation 

As shown in the previous section, larger value of TO3/2 
is allowed compared to the case without f/^j-NLSP. In 
the newly allowed parameter region, the MSSM-LSP 
decays mainly into i>i^, and decays into the grav- 
itino. Since the gravitino is produced with large ve- 
locity dissipation at the late universe, it behaves as 
a warm dark matter, and as a result, may affect the 
structure formation of the universe. In this section, we 
consider the constraints from the structure formation. 



vr decays to ipfj, mainly through two-body process, 
i^R V'ai^a- We find that the lifetime of i>R to^ is 10^- 
10^ seconds for 7773/2 = 0.1-100 GeV with mo,, = 100 
GeV. With the use of t£>^ , we calculate free-streaming 
length Afs = Jr^'^ d,tv{t) /a{t), where v{t) is the ve- 
locity of the i/jfj^, a{t) is the cosmic scale factor, and 
^EQ is the time of the matter-radiation equality. As a 
result, we found Aps — 6 Mpc when m^;^ = 100 GeV 
irrespective of the 7773/2. Thus, it indicates that the 
component of the dark matter (i.e., gravitino) from 
i/R decay acts as a warm dark matter (WDM). In ad- 
dition to vr decay, gravitinos are also produced by 
the thermal scattering at the reheating epoch after in- 
flation. The abundance of the V'/j from the scattering 
process is determined by the reheating temperature 
and the 7773/2 [15]. Since the tpfj, from the scattering is 
non-relativistic at the time of the structure formation, 
it acts as a cold dark matter (CDM). Thus, we have to 
consider the constraints from the structure formation 
of the universe on the WDM+CDM scenario. 

Constraints from the structure formation on the 
WDM-I-CDM scenario are studied in recent works 
[TT] . According to these studies, it turns out that the 
matter power spectrum has a step-like decrease around 
the free-streaming scale of the WDM component, k ^ 
27r/AFS. This fact can be understood intuitively, be- 
cause only the power spectrum of the WDM compo- 
nent dumps at the scale. On the other hand, the power 
spectrum is estimated from the observations of the 
cosmic microwave background T, 3 , the red shift sur- 
veys of galaxies [2j , and so on. Even though the power 
spectrum has been experimentally determined accu- 
rately, an ambiguity still remains when k~^ is around 
1 Mpc. Therefore, it is not clear whether the step-like 
decrease exists or not, if it is small enough. In this arti- 
cle, we put the conservative constraints on our model: 
the power spectrum, to be more precise, the magnitude 
of the step-like decrease, should be within the range of 
the 95% confidence level of the observational data [Tl. 
This condition gives the upper bound on the portion 
of the WDM component. 
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In our model, the energy density of the dark matter 
is composed of two components, puM = P3/2 + 
where ^3^2 P3/2 energy densities of grav- 

itino produced by the decay and by the thermal scat- 
tering processes, respectively. Introducing the fraction 
of WDM component /, we rewrite pdm sls 

PDM - P3/2 + P3/2 

— /PpurcWDM + (1 ^ /)PpurcCDM, (3) 

where PpureWDM and PpureCDM are the energy densities 
of pure WDM and CDM scenario, respectively. Con- 
sidering the adiabatic density fluctuation, we calculate 
the power spectrum. As a result, wc finally get / < 0.4. 
In terms of the density parameter, the constraints in- 
dicates i^3°9 ^ 0.4J?DM, which gives the upper bound 
on the gravitino mass as 7713/2 40 GeV and TO3/2 ^ 4 
GeV for the cases of Cmodci — 1 and 10, respectively 
(the light shaded regions in FigH]). 

The constraints do not depend highly on the detail 
of the observational data. In fact, in other recent ob- 
servations, it is mentioned that the observational error 
on the power spectrum is around 15% plll6j. leading 
to the constraint as / < 0.2, which is of the same order 
of magnitude as the result above. Therefore, our sim- 
ple analysis is expected to give the viable constraint of 
the structure formation on the WDM-I-CDM scenario. 



5 Conclusions 

In this paper, we have studied the cosmological impli- 
cations of the gravitino LSP scenario with the right- 
handed sneutrino NLSP in the framework where neu- 
trino masses are purely Dirac-type. In the case that 
MSSM-LSP is Bino-like ncutralino, it mainly decays 
into the right-handed sneutrino with the lifetime ~ 
10^-10'^ seconds in the wide range of the parameter 
region. Though the MSSM-LSP is long-lived, no visi- 
ble particles are produced in the leading process, thus 
constraints from the BBN scenario is drastically re- 
laxed compared to the case without the right-handed 
sneutrino NLSP. With the quantitative analysis of the 
BBN constraints, we have found the new allowed re- 
gion, 0.1 GeV < 7773/2 < 40 GeV, when m^^ = 100 
GeV. In this region, the BBN constraints give the up- 
per bound on the Bino mass as 777^ < 200-400 GeV, 
which mainly comes from hadronic four-body decays. 
On the other hand, the upper bound on the gravitino 
mass is given by the constraints from the structure for- 
mation of the universe. In our scenario, some part of 
the gravitino is produced by the decay of right-handed 
sneutrino at the late universe. As a result, the grav- 
itino freely streams in the universe and acts as a WDM. 
The gravitino is also produced from thermal scatter- 
ing processes, which acts as a CDM. Taking the CDM 
contribution into account, we have considered the con- 
straints on the WDM-f-CDM scenario from the obser- 
vations of (small scale) structure formation, and finally 
found the upper bound on the gravitino mass. 
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